The Uncorrelated Universe: Statistical Anisotropy and the Vanishing Angular 
Correlation Function in WMAP Years 1-3 
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The large-angle (\ow-£) correlations of the Cosmic Microwave Background (CMB) as reported by 
the Wilkinson Microwave Anisotropy Probe (WMAP) after their first year of observations exhibited 
statistically significant anomalies compared to the predictions of the standard inflationary big-bang 
model. We suggested then that these implied the presence of a solar system foreground, a systematic 
correlated with solar system geometry, or both. We re-examine these anomalies for the data from 
the first three years of WMAP's operation. We show that, despite the identification by the WMAP 
team of a systematic correlated with the equinoxes and the ecliptic, the anomalies in the first- 
year Internal Linear Combination (ILC) map persist in the three-year ILC map, in all-but-one case 
at similar statistical significance. The three-year ILC quadrupole and octopole therefore remain 
inconsistent with statistical isotropy - they are correlated with each other (99.6%C.L.), and there 
are statistically significant correlations with local geometry, especially that of the solar system. The 
angular two-point correlation function at scales > 60 degrees in the regions outside the (kpO) galactic 
cut, where it is most reliably determined, is approximately zero in all wavebands and is even more 
discrepant with the best fit ACDM inflationary model than in the first-year data - 99.97%C.L. for 
the new ILC map. The full-sky ILC map, on the other hand, has a non- vanishing angular two-point 
correlation function, apparently driven by the region inside the cut, but which does not agree better 
with ACDM. The role of the newly identified low-^ systematics is more puzzling than reassuring. 

PACS numbers: 98.80.-k 



I. INTRODUCTION AND RESULTS 

Approximately three years ago, the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) team reported [1-5] 
the results of its preliminary analysis of the satellite's 
first year of operation. While the data is regarded as a 
dramatic confirmation of standard inflationary cosmol- 
ogy, anomalies exist. 

Among the unexpected first-year results, there is a nat- 
ural division into several classes: (1) a lack of large angle 
correlations [1], and violations of statistical isotropy in 
the associated multipoles — 1 = 2 and 3 [6, 7]; plus 
weaker evidence for a violation of statistical isotropy at 
t = 6 and 7 [8]; (2) statistically anomalous values for the 
angular power spectrum, Ce in at least three I bins: a 
trough at I = 22(20 - 24), a peak at I = 40(37 - 44), 
and a trough at I = 210(201 - 220) [3]; (3) hemispheric 
asymmetries in the angular power spectrum over a wide 
range of I [9] ; (4) unexpectedly high cross-correlation be- 
tween temperature and E-mode polarization (TE) at low 
I [1], interpreted by the WMAP team as evidence for a 
large optical depth and thus for very early star formation 
[5]; and a discrepancy between the observed TE angular 
power at the largest scales and the best fit concordance 
model [10]. 

In this paper, we focus on the large angle anoma- 
lies of the microwave background as measured in the 
satellite's first three years of operation, and recently re- 
ported by the WMAP team. In particular, we look at 



how the full three- year WMAP results [11-14] (hereafter 
WMAP 123) are similar to or different from those ob- 
tained from the analysis of first-year data (WMAP1). 
Given the high signal-to-noise of the first-year WMAP 
maps over a large number of well distributed pixels, in the 
absence of newly-identified systematic effects, one would 
not expect significant changes in the low-<? multipoles 
between the one-year and three-year data. The recent 
detailed analysis of [15] quantified this expectation, but 
suggested that the quadrupole, being anomalously low, 
is the least robust of the low multipoles. 

In Section II we repeat our previous analysis of the 
CMB data using the multipole vector framework [7, 16]. 
In Section III we recall the lack of angular correlation 
at large angles, as measured by COBE-DMR [17] and 
confirmed by WMAP1 [5], but not investigated in the re- 
cent analysis of the WMAP team. We therefore include 
our own preliminary analysis of the angular two-point 
correlation function in this paper. We discuss the conse- 
quences of our findings in Section IV. 

For the multipole vector analysis we compute the 
quadrupole and octopole multipole vectors and associ- 
ated area vectors [18] of the WMAP 123 full sky map 
(ILC123) [19], and compare them with those derived 
[7, 16] from three full sky maps based on WMAP1 
data — WMAP's own ILC1 map [20]; Tegmark, de 
Oliveira-Costa and Hamilton's (TOH1) map [21]; and 
the Lagrange multiplier Internal Linear Combinations 
map (LILC1) [22]. We look at the correlations among 
the quadrupole and octopole area vectors, and con- 
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firm that in the ILC123 these are anomalously aligned 
(>99.5%C.L.), as found [7, 16] in the ILC1, TOH1 and 
LILC1 (99.4%C.L.; confirming and strengthening earlier 
results in [6]). We confirm that, as for the first-year maps, 
these alignments are strengthened by proper removal of 
the Doppler contribution to the quadrupole [7]. 

Having established the mutual alignment of the four 
planes defined by the quadrupole and octopole, we re- 
call that in the WMAP1 maps it was found [7, 16] that 
these planes were aligned not just among themselves, but 
also with several physical directions or planes. As in the 
previous work we study these alignments both with and 
without accepting the internal correlations between the 
quadrupole and octopole as given. We find that the align- 
ments generally retain the same significance as in the 
first-year all-sky maps. Only the alignment between the 
ecliptic plane and the quadrupole and octopole planes, es- 
pecially given the latter's internal correlations, is notice- 
ably weakened. Meanwhile, there persist further corre- 
lations of the ecliptic and the quadrupolc-plus-octopolc, 
which rely strongly on the alignment of the quadrupole 
and octopole planes with the ecliptic. 

Finally we also examine the map of the combined 
quadrupole and octopole and observe that, as for 
WMAP1, the ecliptic plane traces a zero of the map 
over approximately one-third of the sky, and separates 
the two strongest cxtrema in the southern ecliptic hemi- 
sphere from the two weakest extrema in the northern 
hemisphere. 

Given the evidence for the lack of statistical isotropy 
on large scales we provide a careful description of the 
definition of the angular two-point correlation function 
means and connections between different possible defi- 
nitions of this function. Our analysis shows that the 
angular two-point correlation function is strikingly de- 
ficient at large angles - in fact, the deficit of power is 
even more significant in the new data than in the first 
year data (99.97%C.L. for the ILC123 in the region out- 
side the kpO galaxy mask). Moreover, we find that the 
quadrupole and octopole power computed directly from 
the cut-sky maps are inconsistent with the quadrupole 
and octopole computed from the maximum likelihood es- 
timators (the latter then being used for all cosmological 
analyses). In the conclusions section we discuss the im- 
plications of these findings. 



II. MULTIPOLE VECTORS 

Several large-angle anomalies or putative anomalies in 
the CMB data have been pointed out and discussed ex- 
tensively in the literature [6, 7, 9, 10, 16, 18, 23-35]. 
Similarly, several novel methods proposed to study the 
statistical isotropy or Gaussianity of the CMB have also 
been introduced and discussed in the literature [36-39]. 
Here we focus on our previous studies [7, 16], namely 
the alignment of the quadrupole and octopole with each 
other and with physical directions on the sky as revealed 



by the multipole vector formalism [18]. 

The £-th multipole of the CMB, Ti, can, instead of be- 
ing expanded in spherical harmonics, be written uniquely 
[18, 40, 41] in terms of a scalar A^ which depends only 
on the total power in this multipole and I unit vectors 
{{j(^>*)|j = 1 ; ...,£}. These "multipole vectors" encode all 
the information about the phase relationships of the ae m . 
Heuristically, 

T t *Awf[(vM.e), (1) 

i=l 

where v^-^ is the i th multipole vector of the £ th mul- 
tipole. (In fact the right hand side contains terms with 
"angular momentum" £—2,^ — 4, ... These are subtracted 
by taking the appropriate traceless symmetric combina- 
tion, as described in [7, 18] and [16].) Note that the signs 
of all the vectors can be absorbed into the sign of 
so one is free to choose the hemisphere of each vector. 
Unless otherwise noted, we will choose the north galactic 
hemisphere when quoting the co-ordinates of the multi- 
pole vectors, but in plots we will show the vector in both 
hemispheres. 

A. Multipole Vectors: WMAP1 vs. WMAP123 

Multipole vectors are best calculated on cleaned full 
sky maps. The ILC1, TOH1, and LILC1 are all 
minimum- variance maps obtained from WMAP's single- 
frequency maps, but differ in the detailed implementation 
of this idea. These full-sky maps may have residual fore- 
ground contamination, probably mainly due to imperfect 
subtraction of the Galactic signal. They also have com- 
plicated noise properties [2] that make them less than 
ideal for cosmological tests. While one can, in principle, 
straightforwardly compute the true (full-sky) multipole 
vectors from the single-frequency maps with a sky cut, a 
cut larger than a few degrees across will introduce signif- 
icant uncertainty in the reconstructed full-sky multipole 
vectors and consequently in any statistics that use them 
(see Section 7 and Figure 7 of [16]). 

The multipole vectors for i = 2 and 3 for the Dopplcr- 
quadrupole-corrected ILC1, TOH1, LILC1 and ILC123 
maps are given in Table I and plotted in Figure 1. 
It is important that the Doppler contribution to the 
quadrupole, inferred from the measured dipole, has been 
removed before computing the multipole vectors. Al- 
though the Doppler-induced piece of the quadrupole is 
a small part of the total power, as shown in [7] and 
[16], failure to properly correct the quadrupole for the 
Doppler contribution results in reduced significance for 
the quadrupole-octopole correlations in all full sky maps. 
Freeman et al. [8] have suggested that both the lack of 
low-£ power and certain reported violations of statistical 
isotropy (in particular any north-south asymmetry) re- 
ported in cut-sky analysis could be due to the use of a 
wrong value for the dipole. However, uncertainty in the 
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TABLE I: Multipole vectors, t)^' 1 - 1 , and area vectors w ( - e ' t ' i \ for the quadrupole and octopole in Galactic coordinates (l,b). 
(Magnitudes are also given for the area vectors.) All vectors are given for the ILC123, ILC1, TOH1, and LILC1 maps after 
correcting for the kinetic quadrupole. 
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dipolc seems not to significantly affect the correlations 
in quadrupole and octopole multipole vectors discussed 
here, as in the full-sky maps the identification of residual 
dipoles can be done uniquely. 

We have found that the area vectors 



;i,3) = fiW.i) x y(i,3)_ 



(2) 



are often more useful than the multipole vectors them- 
selves for statistical comparison. As discussed in [16], 
under certain circumstances relevant to the observed 
quadrupole and octopole, these are closely related to 
the maximum-angular-momcntum-dispersion-axes of [6] . 
They are also related to the Land and Magueijo triad 
[33], at least for the quadrupole, for which that triad is 
uniquely defined. The area vectors for ILC123 are also 
to be found in Table I. 

As expected, the octopole multipole vectors, and con- 
sequently the octopole area vectors, are largely un- 
changed from ILC1 to ILC123. Somewhat unexpectedly 
(although see [15]) the quadrupole vectors have changed. 
This can all be seen quite clearly in Figure 1. To quantify 
the changes in multipole vectors we have computed the 
five dot products between the corresponding year 1 and 
year 123 quadrupole and octopole multipole vectors 
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ILC1 ' ' U ILC123- 



(3) 



Similarly, to quantify the changes in area vectors, we have 
computed the dot products between the old and new area 



vectors 



^(t;i,j) = W ILC1 ' W ILC123 



(4) 



\W 



ILC1 



J ILC123l 



and also the ratios of lengths between year 1 and year 
123 corresponding area vectors 



I W ILC1 I 
l u; ILC123l 



(5) 



These are all found in Table II, together with compari- 
son values for the dot products of vectors in the ILC1 and 
TOH1 maps. In the case of the octopole the change from 
ILC1 to ILC123 is comparable to the differences among 
different full sky maps constructed with the WMAP1 
data but the change from ILC1 to ILC123 is consider- 
ably larger for the quadrupole. 

The change in the ILC quadrupole is partly due to 
WMAP's identification of a systematic effect in the model 
of the radiometer gain as a function of time [11] and 
partly due to the correction of a bias of the ILC map 
[12]. Let us for the convenience of the reader paraphrase 
the details given in [11] and [12] to explain both effects. 

The radiometer gain is the voltage difference measured 
result of changing sky-temperature differences as 
the satellite sweeps the sky. Besides the sky tempera- 
ture, the radiometer gain depends on the temperature 
of the optical and electronic components involved, espe- 
cially the receiver box, which houses the radiometers and 
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FIG. 1: The 1 = 2 (top panel) and I = 3 (bottom panel) multipoles from the ILC123 cleaned map, presented in Galactic 
coordinates, after correcting for the kinetic quadrupole. The solid line is the ecliptic plane and the dashed line is the supergalactic 
plane. The directions of the equinoxes (EQX), dipole due to our motion through the Universe, north and south ecliptic poles 
(NEP and SEP) and north and south supergalactic poles (NSGP and SSGP) are shown. The multipole vectors are plotted as 
the solid red symbols for t — 2 and solid magenta for I — 3 (dark and medium gray in gray scale versions) for each map, ILC1 
(circles), ILC123 (triangles), TOH1 (diamonds) and LILC1 (squares). The open symbols of the same shapes are for the normal 
vectors for each map. The dotted lines are the great circles connecting each pair of multipole vectors for the ILC123 map. For 
1 = 3 (bottom panel), the solid magenta (again medium gray in the gray scale version) star is the direction of the maximum 
angular momentum dispersion axis for the ILC123 octopole. 
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FIG. 2: The quadrupole vectors (filled triangles) and normals (open triangles) of the yrl-yrl23 difference maps for the V band 
(white), W band (black) and the ILC (grey). The background shows the ILC123 quadrupole in the same coordinates as in 
Figs. 1 and 3. 



their electronics. The satellite's temperature and with it 
the temperature of the receiver box changes periodically 
with season and shows a slight long-term increase due to 
the degrading of WMAP's sun shields by the Sun's UV 
radiation [11]. Consequently seasonal modulation and 
longer-term temperature drift of the receiver box must 
be included in a radiometer gain model, which is needed 
to calibrate the sky maps. The gain model is fitted to 
the hourly measurement of the CMB dipole, and conse- 
quently the quality of that fit also depends on the pre- 
cision to which the value and orientation of the CMB 
dipole are known. While in the calibration of WMAP1 
the COBE/DMR dipole has been used, the WMAP123 
analysis relies on the WMAP1 dipole. Having three years 
of data and an improved CMB dipole measurement at 
hand enabled the WMAP team to significantly improve 
the WMAP123 gain model over the WMAP1 data set. 
The effect of the improved gain model would be expected 
to be correlated with the ecliptic plane, the equinoxes 
(seasons) and the CMB dipole, due to the reasons given 
above. The seasonal effect must show up most promi- 
nently in the quadrupole (as there are four seasons). The 
WMAP team estimates that the residual error on the 
quadrupole power is AC2 = 29^K 2 . 

For the synthesis of the ILC (full-sky) map an ad- 
ditional reported systematic correction ("galaxy bias") 
adds to the map a quadrupole and octopole aligned with 



the galaxy to correct for the fact that a minimal vari- 
ance reconstruction of the full sky tends to maximize 
the cancellation between the underlying signal and any 
foregrounds, which are dominantly galactic [12]. The 
WMAP team estimates that the residual error of the full- 
sky ILC after bias correction is < 5{iK at larger than 10 
degrees. This bias correction apparently affects the over- 
all quadrupole power more than the correction to the 
gain model does, but affects the quadrupolar multipolc 
vector structure less. 

In particular, we have checked that it is the new 
gain model that slightly changes the orientation of the 
quadrupole plane (defined by tu^ 2 ' 1,2 )) in the individual 
band maps and in the ILC map. We determined the mul- 
tipole vectors of the lyr-123yr difference full-sky maps 
of the V and W bands. As shown in Fig. 2, we find 
that the quadrupole vectors of the difference maps of 
the bands are quite close to the ecliptic. The associ- 
ated normals are close both to the ecliptic poles, and 
to one of the maxima of the ILC 123 quadrupole within 
the ILC123 quadrupole plane. (Indeed they lie between 
them.) Thus, the lyr-123yr difference maps for the V 
and W bands show a quadrupole that is orthogonal to 
the ILC123 quadrupole and parallel to the ecliptic. Espe- 
cially in the V band (supposed to be the least foreground 
contaminated band), one of the quadrupole vectors of the 
difference map not only lies in the ecliptic plane, but is 
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TABLE II: Comparisons between ILC1 and ILC123 multipole 
vectors and area vectors. Three quantities are tabulated: (a) 
dot products dS l between the multipole vectors from the 
ILC1, t)jLci an d the corresponding vectors from the ILC123, 
*ilci23. for e = 2 through I = 7; (b) dot products A (£;iJ) 
between the normalized area vectors from the ILC1, w^ci 
and the corresponding area vectors from the ILC123, Wj^cimi 
for 1 = 2 and i = 3; (c) ratios r^' 1 '^ between the magnitude of 
the area vectors from the ILC1 \v$£qi | and the corresponding 

magnitudes IwjlcimI fr° m the ILC123,for 1 = 2 and I = 3. 
For comparison, we show in column 3 the same quantity for 
ILC1 versus TOH1. 
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also very close to the equinoxes. 

The additional effect of the ILC bias correction on the 
quadrupole orientation can be estimated by analyzing the 
difference map ILC1-ILC123. Now one of the quadrupole 
vectors is very close to the equinoxes and the quadrupole 
area is close to the quadrupole areas from the V and 
W lyr-123yr difference maps. Compared with the differ- 
ence maps of the individual bands, the second quadrupole 
vector is at higher galactic latitude, which is the signa- 
ture of the admixture of a galactic signal from the ILC 
bias correction (a galactic signal is expected to have its 
quadrupole vectors at high galactic latitudes [16]). 

WMAP's identification of an ecliptic and equinox cor- 
related effect in the radiometer gain is an impressive 
demonstration that the method of multipole vectors is 
very sensitive to systematic effects and/or unexpected 
physics — and has now proven to be a powerful tool for 
modern CMB data analysis. 

On the other hand, we think it is quite disturbing that 
the signature of the systematic correction (the difference 
quadrupole vectors and area) seems to be anti-correlated 
(all the vectors are orthogonal to each other) with the 
supposedly cosmic signal. If both the gain model correc- 
tion and the bias correction would be all that is needed to 
get rid of any non-cosmic signal, one would expect that 
the cosmic signal is uncorrelated with the correction. We 
think that these observations imply the presence of rather 
odd conspiracies between underlying signal and system- 
atics. We suggest that there are additional systematics 
or foregrounds still to be identified. 



B. Multipole Vectors: Violations of Statistical 
Isotropy 

In [7] and [16] we pointed out a number of unusual 
properties of the quadrupole and octopole apparent in 
the WMAP1 full-sky maps (other authors have analyzed 
the WMAP sky using the multipole vectors, e.g. [42-47]. 
Here we review these properties for the ILC123 maps 
as illustrated in Figures 1 and 3. We refer the reader 
to Section 4 of [16] for further details concerning ILC1, 
TOH1 and LILC1. 

1. The Queerness of the Quadrupole (cf. Fig. 1 top 
panel). The changes in the quadrupole multipole 
vectors are noticeable. 

(a) The normalized quadrupole area vector 
^(2 ; i,2) gt -ij Y[ es near the ecliptic plane but has 
moved slightly away from the plane. 

(b) t/)( 2;1 > 2 ) appears somewhat aligned with both 
the dipole and the equinoxes but has also 
moved away from these locations. 



2. The Oddness of the Octopole (cf. Fig. 1 bottom 
panel). The octopole has not changed beyond ex- 
pectations. 
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FIG. 3: The I — 2 + 3 multipoles from the ILC123 cleaned map, presented in Galactic coordinates. This is a combination of 
the two panels of Fig. 1 with only the multipole vectors for the ILC123 map shown for clarity. The solid line is the ecliptic 
plane and the dashed line is the supergalactic plane. The directions of the equinoxes (EQX), dipole due to our motion through 
the Universe, north and south ecliptic poles (NEP and SEP) and north and south supergalactic poles (NSGP and SSGP) are 
shown. The t = 2 multipole vectors are plotted as the solid red (dark gray in gray scale version) triangle and their normal 
is the open red (dark gray in the gray scale version) triangle. The I = 3 multipole vectors are the solid magenta (medium 
gray in gray scale version) triangles and their three normals are the open magenta (medium gray in the gray scale version) 
triangles. The dotted lines are the great circles connecting the multipole vectors for this map (one for the quadrupole vectors 
and three for the octopole vectors). The solid magenta (again medium gray in the gray scale version) star is the direction to 
the maximum angular momentum dispersion for the octopole. 



(a) Two of the normalized octopole area vectors, 
$(3:1,2) an( j ^(3;3,i)^ ij e near ecliptic plane. 

(b) The remaining normalized octopole area vec- 
tor, uj( 3 ' 2 ' 3 \ lies very near the supergalactic 
plane, and about 10° from the Galactic pole. 

(c) The octopole appears somewhat planar — the 
three normalized octopole area vectors cluster 
somewhat more than expected. 

3. The Remarkable Relation of the Quadrupole and 
Octopole (cf. Fig. 3). For the most part the re- 
lationship of the quadrupole and octopole hasn't 
changed. 

(a) The normalized quadrupole area vector is 
aligned with the normalized octopole area vec- 
tors in the sense that it lies "in the middle" 
of these vectors. Furthermore, it now lies very 
close to the octopole's maximum angular mo- 
mentum dispersion direction [6] . 

(b) The ecliptic carefully traces a zero of the com- 
bined map, almost perfectly over the entire 



hemisphere containing the Galactic center, 
and relatively well over the antipodal hemi- 
sphere. 



(c) Two of the three extrema south of the eclip- 
tic are clearly stronger than any north of the 
ecliptic. The strongest of the three northern 
extremum is only comparable in strength to 
the weakest southern extremum. 



The above observations are purely qualitative. Their 
statistical significance must be evaluated. To this end we 
developed several statistical tests for application to the 
first-year maps; we now repeat these tests for ILC123. 
Although we have also applied some of these tests to 
higher multipoles (with possibly interesting results), we 
defer consideration of higher multipoles to a future paper. 
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TABLE III: Comparison of the values of the A; (the dot prod- 
ucts of the quadrupole area vector with each of the three oc- 
topole area vectors). Ai and S(Ai) = (Ai + A2 + As)/3 are 
tabulated for ILCl, TOHl, LILCl and ILC123, as well as the 
percentile rank (out of 100) of S(Ai) among a suite of 10 
Gaussian random statistically isotropic skies. For the ILC123 
map we show the results with and without the DQ-correction. 





ILCl TOHl LILCl ILC123 ILC123 








uncorr 




A l 


0.9251 0.8509 


0.7798 


0.716 


0.7666 


A 2 


0.7864 0.7829 


0.7434 


0.757 


0.7924 


A 3 


0.6077 0.7616 


0.7229 


0.626 


0.7221 


S(Ai) 


0.7731 0.7985 


0.7487 


0.700 


0.7604 


-P[S(Ai)] 


0.289 0.117 


0.602 


1.831 


0.433 



1. Correlations Among Multipole Vectors 

We first consider the alignment between the 
quadrupole area vector and the three octopole area vec- 
tors. This alignment can be measured by the magni- 
tudes of the dot products between w( 2;1 < 2 ) and each of 
yj(3;»,j) _ w e define these, according to the ordering from 
the largest to the smallest for the first-year full-sky maps 
(ILCl, TOHl, LILCl): 





1^1,2) 


■ w { 


A 2 = 


N (2;i,2) 


■ w { 


A 3 = 


|^(2;i,2) 


■ w l 



A natural choice of statistic which defines an ordering 
relation on the three dot products is [7, 16, 43] 

S(A i )= 1 -f2A i . (7) 

In Table III the values of the Ai and S(Ai) are shown for 
ILCl, TOHl, LILCl and ILC123. As expected, the dif- 
ferences between ILCl and ILC123 appear comparable to 
the systematic errors in making full sky maps, as mea- 
sured by the differences between the ILCl, TOHl and 
LILCl values. However, in contrast to the first year maps 
the ILC123 obeys the ordering A 2 > A\ > A3, which is 
due to the fact that the area vector of the quadrupole 
has changed significantly. 

We have compared the values of the A, against 10 5 
Monte Carlos of Gaussian random statistically isotropic 
skies with pixel noise (as in [18]). Of course, WMAP1 
data is compared to Monte Carlo simulations with the 
one-year pixel noise, while for WMAP123 data we use the 
appropriate three-year pixel noise. The percentile rank 
V of the S(Ai) values for each of the four full sky maps 
among the 10 5 associated Monte Carlos is also listed in 
Table III. The clear inference to be drawn is that the 
quadrupole and octopole show a statistically significant 
correlation with each other in the ILC123 map, just as 



they do in the WMAP first-year all-sky maps. This cor- 
relation is at the 99.6%C.L. for ILC123, within the range 
seen for the WMAP1 maps (99.4 - 99.9%C.L). 

Also included in Table III arc the values of the A; 
and the associated statistics for an ILC123 from which 
the Doppler contribution to the quadrupole has not been 
removed (ILC123 uncorr). As for WMAP 123, we see 
that the correlation between quadrupole and octopole is 
significantly stronger in the properly corrected map than 
in the uncorrected map. 

2. Correlations of area vectors with physical directions 

So far we have looked only at the internal correla- 
tions between the CMB quadrupole and octopole. We 
have previously [7, 16] examined the correlation of the 
quadrupole and octopole area vectors with various phys- 
ical directions, and inferred that the quadrupole and oc- 
topole are also correlated not just with each other but 
also with the geometry of the solar system. We now re- 
examine these correlations quantitatively for WMAP 123. 

Our measure of alignment of the area vectors with a 
physical direction d is 

S^\d) = ±(\d-w^\ + \d-w^\ (8) 

+ |d. w j(3;2,3)| + | d '. uJ (3;3,l)|y 

In Table IV, we display the probabilities of the S^ 4 ' 4 \d) 
statistic for the various WMAP1 and WMAP123 full-sky 
maps computed for: the ecliptic plane, the galactic poles 
(NGP), the supergalactic plane, the CMB dipole and the 
equinoxes. These directions and planes are the ones most 
obviously connected either to possible sources of system- 
atics (ecliptic plane, CMB dipole and equinoxes) or to 
possible unanticipated foregrounds (ecliptic plane, galac- 
tic plane/pole, supergalactic plane). For the galactic 
poles, the dipole and the equinoxes, these probabilities 
are the percentage of statistically isotropic skies (among 
our usual suite of 10 5 Monte Carlo realizations of sta- 
tistically isotropic skies) that exhibit a higher value of 
gtM) f or that direction; for the ecliptic and supergalactic 
planes, they are the percentage of statistically isotropic 
skies (among the 10 5 realizations) that exhibit a lower 
value of S^ 4 ' 4 ) for the axis of each plane. 

From Table IV, we see that the correlations of the full 
sky map to the dipole, equinox and galactic poles remains 
essentially unchanged in the WMAP 123 at 99.7%C.L., 
99.8%C.L. and 99%C.L. respectively. Similarly the cor- 
relation to the supergalactic plane remains insignificant 
at 85%C.L. The correlation with the ecliptic has declined 
somewhat from 98%C.L. to 96%C.L. 

In Section II B 1 we established the correlation of the 
quadrupole and octopole area vectors with each other. 
One might therefore ask whether the observed correla- 
tions to physical directions could be mere accidents of 
the internal quadrupole-octopole correlations. In [16], 
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TABLE IV: The percentiles of 5 (4,4) (d) from 10 s MC 
maps test for the alignment of known directions with the 
quadrupole and octopole area vectors. The ILC1, TOH1, 
LILC1 and ILC123 maps are studied. For the ILC123 map 
we show the results with and without DQ-correction. 





ILC1 


TOH1 


LILC1 


ILC123 ILC123 










uncorr 




ecliptic 


2.01 


1.43 


1.48 


4.88 


4.11 


NGP 


0.51 


0.73 


0.94 


1.08 


0.89 


SG plane 


8.9 


14.4 


13.4 


10.8 


15.3 


dipole 


0.110 


0.045 


0.214 


0.489 


0.269 


equinox 


0.055 


0.031 


0.167 


0.328 


0.194 



we showed that for the first-year maps, even given the 
"shape" of the quadrupole and octopole their multipole 
structure and mutual orientation - the additional corre- 
lations of their area vectors to the ecliptic plane were 
very significant. We also found that the additional cor- 
relations to the dipole and equinoxes was less significant, 
and the additional correlations to the galaxy and to the 
supergalactic plane were not statistically significant. We 
next revisit the issue for the ILC123. 

To address this question, we compute the probability, 
given the shape of the quadrupole and octopole and their 
mutual alignment, that they would align to the extent 
they do with each direction or plane. Put another way, 
we find the fraction of directions/planes that are bet- 
ter aligned with the observed quadrupole and octopole 
than each of the directions and planes in question. Thus, 
we hold the quadrupole and octopole area vectors fixed 
and compute S^ 4 ' 4 '(d) both for the physical directions 
that may be of interest — the ecliptic poles, the Galac- 
tic poles, the supergalactic poles, the cosmological dipole 
and the equinoxes — and for a large number of random 
directions. We denote this test by S^ 4 ' 4 ) (d\{w}), mean- 
ing that the area vectors {w} are fixed to be the observed 
ones. 

The results are shown graphically in Figure 4, which 
is a histogram of the values of S^ 4,4 \d\{w}) for the ran- 
domly chosen directions, with the values of S^ 4 ' 4 ^ (d\{w}) 
for the particular physical directions shown. The inter- 
esting shape of the histogram, including the spike, is dic- 
tated by the geometry of the particular mutual orienta- 
tion of the quadrupole and octopole. Directions lying 
within the relatively small triangle bounded by the three 
octopole area vectors - which contains the quadrupole 
area vector - contribute to the spike. The existence of the 
spike is therefore a consequence of quadrupole-octopole 
alignment. 

It is clear from Figure 4 that the additional correla- 
tions with the galaxy and with the supergalactic plane 
are not particularly significant. To determine the precise 
significance we compute for each physical direction, the 
fraction of random directions that yield a larger value of 
S( 4 < 4 )(d| {«;}). These are given in Table V. 
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FIG. 4: Histogram of the S^- 4 ' 4 ^ (d\{w}) statistics applied to 
the ILC123 map quadrupole and octopole area vectors and 
a fixed direction or plane on the sky, compared to 10 5 ran- 
dom directions. Vertical lines show the 5* statistics of the 
actual area vectors applied to the ecliptic plane, NGP, super- 
galactic plane, dipole and equinox directions (Table V shows 
the actual product percentile ranks among the random ro- 
tations for ILC1, ILC123 and TOH1). This Figure and Ta- 
ble V show that, given the relative location of the quadrupole- 
octopole area vectors (i.e. their mutual alignment), the dipole 
and equinox alignments remain unlikely at about 95% C.L., 
whereas the ecliptic alignment, significant at 98.3-99.8% C.L. 
in year 1, is only significant at the 90% C.L. in year 123. The 
galactic plane and supergalactic plane alignments remain not 
significant. 



TABLE V: The percentiles of S (4A) {d\{w}) for the ecliptic 
plane, NGP, supergalactic plane, dipole and equinox axes 
among a comparison group of 10 random directions (or 
equivalently random planes) on the sky. The ranks are given 
for the ILC1, TOH1, LILC1 and ILC123 maps. 



ILC1 TOH1 LILC1 ILC123 



ecliptic 


1.7 


1.0 


0.2 


10.3 


NGP 


90 


87 


88 


88 


SG plane 


25 


34 


33 


32 


dipole 


94.5 


95.6 


93.8 


93.0 


equinox 


96.4 


96.1 


94.4 


94.0 



In Table V we compare the results obtained from 
ILC123 to those obtained from ILC1, TOH1 and LILC1. 
As expected, the dipolar and equinoctic values for 
S 1 * 4 ' 4 ^ (d\{w}) remain essentially unchanged, in the 93th 
percentile and 94th percentile respectively among all di- 
rections. However, one should probably interpret this 
as just a 86% C.L. and 88%C.L. case respectively for 
additional correlation, because we would have found it 
equally surprising if they had been in the 7th or 6th per- 
centile. The 88th percentile rank of the Galactic polar 



10 



axis similarly offers only a 76% C.L. case for additional 
correlation, while the supergalactic plane is clearly not 
additionally correlated by this measure. 

Our first real surprise is that the rank S^ 4 ' 4 ) (d|{w}) for 
the (normal to the) ecliptic plane has increased from the 
0.2 — 1.7 percentile range of the three WMAP1 maps to 
10 percentile for the ILC123. This is because, although 
the altered quadrupole caused only a small increase in 
the actual value of S^ 4 ' 4 \d\{w}) for the dipole, the value 
moved from below the peak at low S^ 4 ' 4 ) seen in Figure 
4, into the main body of the peak. This still represents 
90% C.L. evidence of additional correlation of the eclip- 
tic, not 80% C.L. evidence. That is because, as we shall 
see in subsection II B 3, the ILC123 possesses other fea- 
tures strongly correlated with the ecliptic, not captured 
by S^ 4 ' 4 \d\{w}) (or even by the area vectors). These 
features would be impossible if the ecliptic S^ 4 ' 4 \d\{w}) 
were in the top 10 percentile of values. 

On the basis of the S^ 4 ' 4 ^(d\{w}) statistic alone, while 
the additional correlation with the ecliptic is the most sig- 
nificant among the physical directions and planes which 
we have considered, it is now of only 90%C.L. signifi- 
cance. However, as we have suggested above, and ex- 
plore below, the S^ 4 ' 4 \d\{w}) statistic captures only a 
part of the strange connection between the quadrupole 
and octopole and the ecliptic. 

Some measure of caution should also be retained in 
dismissing any Galactic, supergalactic, dipole or equinoc- 
tic correlations just because the additional correlations 
are not significant. We do not know whether the cor- 
relations with those directions shown in Table IV are 
"accidental" consequences of the internal correlation of 
the quadrupole and octopole, or whether the correlation 
among the quadrupole and octopole area vectors is an 
"accidental" consequence of their correlation to a physi- 
cal direction. 

Finally, we caution the reader in advance that in Sec- 
tion III we will call into question some of the "improve- 
ments" in the ILC123 over the ILC1 that have most 
likely led to the unexpected decline in statistical signif- 
icance of the additional correlation of the ecliptic with 
the quadrupole and octopole. 

3. Correlations with physical directions not captured by 
area vectors 

The statistics we have studied above do not use all the 
information in a multipolc. The area vectors alone do 
not contain all this information; this is obvious for the 
quadrupole, and true for higher t as well. Furthermore, 
we have considered only the dot products of the area 
vectors with the physical directions and planes, this again 
reduces the portion of information retained. 

In Figure 3, we see very clearly that the ecliptic plane 
traces out a locus of zero of the combined quadrupole and 
octopole over a broad swath of the sky — neatly sepa- 
rating a hot spot in the northern sky from a cold spot 



in the south. Indeed, it seems to separate three strong 
southern extrema from three weaker northern ones. This 
information is not contained in the 5 ,( - 4 ' 4 ^ statistic since it 
depends precisely on an extra rotational degree of free- 
dom about an axis close to the ecliptic plane. As dis- 
cussed at some length in [16], given the observed internal 
correlations of the quadrupole and octopole area vectors 
and given the correlation of those area vectors with the 
ecliptic, the probability of the ecliptic threading along a 
zero curve in the way it does is difficult to quantify ex- 
actly, but is certainly less than about 5%. Furthermore, 
this correlation could not occur if the area vectors were 
pointing toward the ecliptic poles instead of perpendic- 
ular to them. (Hence our assertion that the probability 
of the observed area vectors' correlation with the eclip- 
tic is 10%, not twice that.) This remains as true for the 
ILC123 maps as it was for the WMAP1 full sky maps. 

4- Summary 

To conclude our re-analysis of the quadrupole and oc- 
topole multipole vectors, there is still strong evidence 
for the quadrupole-octopole alignment and/or the align- 
ment of the quadrupole and octopole with the eclip- 
tic/equinox and dipole. The data do not support the sta- 
tistical isotropy of the sky at the largest angular scales. 

III. ANGULAR CORRELATION FUNCTION 

The near vanishing of the two-point angular correla- 
tion function at angular separations greater than about 
60 degrees has been the longest established among the 
anomalies in the CMB [17]. Here we compare the results 
from the one-year WMAP data and the three-year data. 
Since WMAP has not provided an analysis of the angular 
correlation function for the year 123 data we provide our 
own. 



A. What is the "angular correlation function"? 

There are several distinct quantities that are com- 
monly referred to as the "two point angular correlation 
function" and the related "angular power spectrum" . 
Here we clarify how the various definitions depend on the 
assumptions about statistical isotropy and sky coverage 
of the map. 

We have first 

C(ei,e 2 ) = <T(ei)T(e 2 )> (9) 

where (•) represents an ensemble average over the tem- 
peratures T(e.\) and T{e.2) in the particular directions 
e-y and e 2 . (The precise nature of the ensemble when 
we have only one universe to observe is philosophically 
troubling to many, but will be ignored.) If the sky were 
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statistically isotropic, then C(ei, £2) depends only on the 
angle 9 between e.\ and e 2 . This allows us to define our 
first two-point angular correlation function 

C{6) = (T(e!)T(e 2 )) fl) (10) 

where now represents an ensemble average over the 
temperatures T(e{) and T(e.2) in all pairs of directions 
ei and e 2 separated by the angle 9. The average (•)# 
may stand for the full sky or just a finite region of the 
sky. In the case of statistical isotropy any region of the 
sky provides the same answer, however, in the case of a 
statistically anisotropic sky, different regions of the sky 
give rise to different values of C(9). 

We can represent C(6), as a sum involving the Legendre 
polynomials: 

00 

C{6) = — X>^+l)<WMcos0), (11) 
e=o 

C e = 2tt J d{cosd)C(6)P e (cos6). (12) 

These Ce are also our first incarnation of the "angular 
power spectrum" . 

The two point angular correlation function of equation 
(10) is defined as an ensemble average. By definition, we 
can get an unbiased estimate C(9) of C{9) by replacing 
the ensemble average (for the observed part of the sky) 
with an average over the observed sky, even if this is not 
the full sky: 

C(9)^T(e 1 )T(e 2 ) e . (13) 

This average is over all pairs of directions e\ and £2 sep- 
arated by the angle 9 on the actual sky. The more, inde- 
pendent pixels separated by an angle 9 measured on the 
sky, the more accurate the estimate; but even if the sky 
is cut, the estimate is unbiased. C{9) can also be written 
as a sum over Legendre polynomials: 

00 

C(9) = — ^(^+1)0^(008 9). (14) 

e=o 

An alternative to performing the average over the ob- 
served sky and then a Legendre polynomial expansion to 
obtain the Ce is to do a maximum likelihood estimate 
(MLE) to obtain C| 1LE . This has several advantages. It 
allows one to account for noise; for the prior probability 
distribution of the Ce (in a statistically isotropic model 
for the origin of fluctuations, such as inflation); and for 
residual galactic uncertainty [48] (by marginalizing over 
foreground templates). Because the MLE method gives 
a distribution of likelihoods as a function of the Cj vlLE , 
one has a measure not just of the most likely value of 
each Ce but also of the uncertainty in that value. 

Some caution is necessary in using the MLE. MLE as- 
sumes that each T(e^) is Gaussian distributed. It also 
assumes statistically isotropic noise. For I < 10, the 



WMAP team quotes [12] the peak in the likelihood dis- 
tribution of C^ LE as the value of the "angular power 
spectrum" for that I. 

There is another closely related quantity, which is of- 
ten called the two-point angular correlation function, but 
which is not identical to the ensemble average (10) except 
in certain circumstances. Consider the usual spherical 
harmonic expansion of the CMB sky 

00 00 1 

T(9,cf>) = Y / Te = Y, E a *™ Y e™ & ■ ( 15 ) 
e=o e=o m=-i 

The ensemble average of the product of ae m , 

Cemt'rn' = (a*t m a £' m') , (16) 

encodes the same information as C(ei,e2). The sky is 
statistically isotropic if and only if 

Cemt'm' — Ce5u<S mm '. (17) 

In this case one can show that the angular power spec- 
trum Ce are identical to the coefficients in equation (11), 
i.e. 

statistical isotropy => (p t = Ce) . (18) 
Similarly, one can define a function based on the Ce, 

1 oo 

C(ff) = —J2(^+l)CePe(cos9). (19) 
n 1=0 

Unfortunately, we cannot compute the Ce from what 
we measure since the Ce are properties of the underlying 
ensemble of which the sky is but one realization. What 
we can compute are estimators of Ce- For example, given 
the full sky, we could measure the ae m from the full sky 
and then define 

1 1 

C ^WTT £ |a ^' 2 - (20) 

m— — e 

Obviously this definition is motivated by statistical 
isotropy, but is not restricted to that case. In the case of 
statistical isotropy it is an unbiased estimator of Ce- Ad- 
ditionally, if Gaussianity holds, this is the best estimator 
with cosmic variance Var(C^) = 2C|/(2£ + 1). If Gaus- 
sianity is violated the cosmic variance may be larger. 

The estimator (20) can be easily adapted to cut skies. 
Other similar estimators exist (and have improved per- 
formance, for example in a noisy environment [12]). If I 
is large, Ce will be a good estimate of the angular power 
spectrum Ce- These Ce, however precisely obtained, are 
often called the pseudo-C^. 

Finally, one can define the function based on the Ce, 

-. oo 

C(9) = —Y,W+±)CePe(cos9). (21) 
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TABLE VI: The four classes of pairs of quantities often called 
the "two-point angular correlation function" and the "angu- 
lar power spectrum". Angle brackets represent ensemble av- 
erages, an overbar represents an average over the sky or some 
portion thereof. LPT stands for Legendre polynomial trans- 
form, LPS for Legendre polynomial series expansion. Tj is 
the temperature in the ith sky pixel, at m is the coefficient 
of Y( m in a spherical harmonic expansion of the temperature 
field on the sky. In the case of Ce, the definition actually 
represents the canonical, but not necessarily best, estimator 
of Ce, calculated from the spherical harmonic coefficients of 
the temperature field. 
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In the case of statistical isotropy and Gaussianity the cos- 
mic variance of the correlation function estimated from 
full skies becomes 

oc 

Var(C(0)) = £ (2 ^+ l)CfP, 2 (cos0). (22) 

1=0 

If the Ce and C(9) are inferred from a full sky measure- 
ment, then one can show that (up to pixelization errors) 

full sky (Ce = C e ) , (23) 

and equivalently 

full sky => (C(9) = C(9)) . (24) 

The relationships among these quantities are summarized 
in Table VI. 

While the temptation not to differentiate between 
these quantities is understandable, we would argue for 
a preference to reserve for C(9) (and C{9)) the name an- 
gular correlation function (and its estimator) ; and to call 
Ce (and Ce) the angular power spectrum (and its estima- 
tor). 

In the previous Section we presented evidence that the 
low-^ sky is not statistically isotropic. Moreover the ae m 
(or the Ce) are generically not computed over a full sky. 
Therefore the distinctions between C and C , and between 
C and C must be kept in mind when interpreting CMB 
results. In particular, the measurable quantities — C(6) 
and Ce — are not interchangeable, but rather encode, at 
least in principle somewhat different information about 
the map. This means that one must check that any anal- 
ysis procedure has reasonable consequences not just for 
the Ce, but also for C(9). 

Indeed, in the absence of statistical isotropy, many of 
the quantities to which we are accustomed are of uncer- 
tain meaning and dubious value. For example, neither 



C(9) nor Ce characterizes the statistical properties of the 
ensemble of which the CMB is a realization. They are 
averages of quantities that may (partly) characterize the 
ensemble — C(e\, e-i) or Ce-mi'm' — but, in the absence of 
SI, these more complicated quantities cannot be mean- 
ingfully estimated with only one sky. Estimators like 
C(9) and Ce remain of interest primarily to the extent 
that either the sky is almost statistically isotropic (so 
that C(e±, ei) and Ce-mi'm' can be expanded around C(9) 
and CeSee'Smm') or the estimators (C(9) or Ce) return 
unexpected values. 

The question of approximate statistical isotropy may 
be worth exploring. What is certainly true, as we show 
below, is that \C(9)\ remains unexpectedly small at large 
angles in WMAP123, just as it was for COBE-DMR 
and WMAP1. One must therefore be extremely cautious 
about any reconstruction method that regards the Ce as 
fundamental, independent variables, while ignoring the 
properties of C(9). 

To summarize, while CMB cosmologists have a long- 
standing interest in and preference for the angular power 
spectrum Ce as the relevant variables for statistical anal- 
ysis, given the absence of statistical isotropy for low I 
these variables are at best flawed. Therefore one must 
be careful to distinguish between Ce and Ce and their es- 
timators; one must be cautious that Ce and Ce may not 
characterize the observed microwave maps; and one must 
pay careful attention to the properties of C(9). 

B. The measured C(0) 

We can readily explore C(9), the angular correlation 
function, since there is indeed a very large number of in- 
dependently measured pixels on the WMAP sky. C(9) 
was first measured (at large angles) using the Cosmic 
Background Explorer's Differential Microwave Radiome- 
ter (COBE-DMR) [17], and found to be anomalously 
small in magnitude at large angles in all COBE bands. 

The WMAP team [1, 5] using their first year of data 
presented the "angular correlation function", but this 
was obtained from the Ce (and so, was actually C(9), 
not C(9)). WMAP did not present any angular corre- 
lation function in their third-year papers, although, as 
mentioned above, their "angular power spectrum" at low 
I were actually maximum likelihood estimates of Ce, i.e. 
the coefficients of a Legendre polynomial series expansion 
of a MLE of C(9). 

In Figure 5 we show C(9) = T(ei)T(e2) e (computed 
in pixel space) for the Q, V and W band maps masked 
with the kpO mask. (See [1] for a discussion of the vari- 
ous galactic masks.) The K and Ka bands have a large 
and obvious galactic contamination even after applying 
the kpO mask, and we do not use them. Also shown 
are the correlation function for the ILC map with and 
without the mask, and the Legendre polynomial series 
of the quoted maximum likelihood estimates of Ce, i.e. 
C MLE (9). Finally, we show C(9) = C{9) for the best fit 
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FIG. 5: Two point angular correlation function, C(ff) = T(ei)T(e2) g computed in pixel space, for three different bands masked 
with the kpO mask. Also shown is the correlation function for the ILC map with and without the mask, and the value expected 
for a statistically isotropic sky with best-fit ACDM cosmology together with 68% error bars. Left panel: year 1 results. Right 
panel: year 123 results. Even by eye, it is apparent that masked year 123 maps have C(6) that is consistent with zero at 9 > 60 
deg, even more so than in year 1 maps. We also show the C(6) computed from the "official" published Ce, which (at £ < 10) 
are the pseudo-ft in year 1, the and MLE Ce in year 123. Clearly, the MLE-based Ce, as well as C(8) computed from the 
full-sky ILC maps, are in significant disagreement with the angular correlation function computed from cut-sky maps. 



ACDM model of the three-year WMAP, and a blue band 
around that showing the range of C(9) that one would 
extract from a full sky, given cosmic variance. The left 
panel shows the year 1 results, while the right panel shows 
year 123 results. 

First of all, we see that the different bands and the 
cut-sky ILC map are all in excellent agreement with each 
other for year 1 and year 123 separately, perhaps with an 
exception of the year 1 Q band map. In each case the 
magnitude of C(9) is very small above 60 degrees and 
below 170 degrees; above 170 degrees there is a peculiar 
slight anti-correlation while the model predicts a large 
positive correlation. Moreover, even by eye it is appar- 
ent that year 123 maps have C(9) that is nearly vanishing 
at 60 < 9 < 170 deg even more precisely than in year 1 
maps. Finally, the full-sky C(9) from the ILC maps is in 
agreement with i.e. C MLE ((9), but in significant disagree- 
ment with the cut-sky value. 

We now compare the smallness of C{9) to that expected 
in simulated maps. We use the statistic proposed by 
WMAP [5] in their first-year papers that quantifies the 
lack of power at 9 > 60 degrees 

r 1/2 2 

S 1/2 = J [C(9)\ 2 d(cos0). (25) 

The WMAP team computed 6*1/2 for the original ILC1 
and found that only 0.15% of the elements in their 
Markov chain of ACDM model CMB skies had lower val- 
ues of S1/2 than the true sky. In addition we define the 



statistic 

■Sfuii = J [C(0)f d(cos0), (26) 

which quantifies the amount of power squared on all an- 
gular scales. 

In Table VII we give the value of Si/ 2 for the one- year 
and three-year band maps, for ILC1 and ILC123, and 
for the best fit ACDM model. The values of S1/2 and 
S'fuii are very low relative to what the ACDM isotropic 
cosmology predicts and are noticeably lower in the year 
123 maps than in the year 1 maps. 

To quantify the probabilities of our statistics, we per- 
form comparison with 20,000 Gaussian random, isotropic 
realizations of the skies consistent with models favored 
by the WMAP data. The angular power spectrum from 
each simulated sky is generated as follows: 

1. pick a model (in order of the highest weight) from 
WMAP's publicly available Markov chains; 

2. generate the angular power spectrum in multipolc 
space corresponding to that model by running the 
CAMB code [49]; 

3. generate a map in Healpix [50] consistent with the 
model; 

4. apply the kpO mask; 

5. degrade to NSIDE=32 resolution, then reapply the 
kpO mask also degraded to NSIDE=32; and 
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TABLE VII: Values of S 1 / 2 and Sf u ii statistics, and their ranks (in percent) relative to Monte Carlo realizations of isotropic 
skies consistent with best-fit ACDM models from WMAP 3 year constraints. All of the maps were subjected to the kpO mask 
and compared to 20,000 Monte Carlo realizations with the kpO mask also applied. The third and second-to-last row are without 
the mask and compared to 1000 full-sky simulations, while the last row shows the ACDM expectation. The boldfaced entries 
show the extremely small values of Si/ 2 statistic for the 3-year cut-sky maps, while the single italic entry shows that the 
full-sky ILC123 map does not have an unusually small C(8). We also show the quadrupole, octopole and hexadecapole power 
Ce estimated from the corresponding C(8); note again the difference between the cut-sky and full-sky values. 
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6. compute C(6) directly from the map, using the un- 
masked pixels. 

In addition, and using the same procedure but with- 
out masking, we compute C{9) from 1000 simulated full- 
sky maps for comparison to full-sky ILC maps. (Since 
P(5i/ 2 ) and -P(<Sfuii) f° r the full-sky maps are not as low 
as for the cut-sky maps, we did not need to go through 
the resource-consuming task of generating 20 times more 
simulations in order to get good statistics.) We have 
checked that applying the inverse-noise weighting to the 
pixels when degrading the map resolution (O. Dore, pri- 
vate communication) leads to negligible differences in the 
recovered C(0) (< lfiK 2 on large scales). Finally we have 
checked that, for the cut-sky case, applying the kp2 mask 
instead of kpO gives similar results. 

Results from the comparisons are given in the third 
and fifth column of Table VII. They indicate that only 
0.04%-0.14% of isotropic ACDM skies (for the Q, V and 
W band maps) give S1/2 that is smaller than that from 
WMAP123. Similarly, the masked ILC123 map has S 1/2 
low at the 0.03% level. Interestingly, even the Sf u u statis- 
tic that includes power on all scales, applied to cut-sky 
maps, is low at about a 1% level. We therefore con- 
clude that the absence of large angle correlations at scales 
greater than 60 degrees is, at > 99.85 %C.L., even more 
significant in year 123 than in year 1. 

We also show the quadrupole, octopole and hexade- 
capole power inferred from C{9) of various maps, as 
well as the ACDM expected values. It is clear that the 
quadrupole and octopole computed from cut-sky C{6) are 
a factor of two to five smaller than those from the full-sky 
maps (and of course the latter are another factor of two 
to four smaller than the most likely ACDM values). This 
result has been hinted at by the WMAP team [12] who 



pointed out that the MLE estimates of the quadrupole 
and octopole (as well as I = 7) they adopted in year 123 
are considerably larger than the pseudo-CV values they 
had adopted in year 1. Here we find that the pseudo- 
Ce estimates are in a much better agreement with the 
cut-sky C(8) than the MLE values; sec Figure 5. 

The current situation seems to us very disturbing. 
Standard procedures for extracting the "angular power 
spectrum" of the CMB uses pseudo-C^. The recent 
WMAP analysis included a maximum-likelihood analysis 
at low I (described above) to correct for the effects of a 
cut sky and to mitigate any remaining galaxy contami- 
nation. But both the connection between pseudo-CV and 
any properties of the underlying ensemble (whatever that 
may mean), and the validity of the maximum- likelihood 
analysis rely on the assumption that the underlying sky 
(i.e. the ensemble) is statistically isotropic. However, 
we (and others) have already provided strong evidence 
that the underlying sky, or at least the all-sky maps, are 
not statistically isotropic. In the absence of statistical 
isotropy none of the computable quantities have any clear 
connection to the statistical properties of the underlying 
ensemble. Indeed, in the absence of statistical isotropy 
each ag m could have its own distinct distribution. Since 
we have precisely one sample from each such distribution 
- the value of ag m on the actual sky (and even that only 
for a measured full sky) - it is impossible to in any way 
characterize these distributions from the observation of 
just one universe. 

We are faced with two choices: cither the universe 
is not statistically isotropic on large scales (perhaps be- 
cause it has non-trivial topology) , or the observed viola- 
tion of statistical isotropy is not cosmological, but rather 
due either to systematic errors, or to unanticipated fore- 
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grounds. In this latter case, it is imperative to get to the 
root of the problems at large angles. 

The disagreement between the full-sky and cut-sky an- 
gular correlation functions, evident in Figure 5 and Ta- 
ble VII, is both marked and troubling. The third and 
second-to-last columns of Table VII show that both the 
quadrupole and the octopole power estimated from the 
cut-sky maps are much smaller than those estimated from 
reconstructed full-sky ILC maps - for each multipole, the 
former is about one half as large as the latter. However, 
it is the the full-sky values which correspond to the esti- 
mated values used in cosmological parameter analyzes. 

It is of utmost concern to us that one begins with in- 
dividual band maps that all show very little angular cor- 
relation at large angles (i.e. \C{6)\ is small), and, as a 
result of the analysis procedure, one is led to conclude 
that the underlying cosmological correlations are much 
larger. What is even more disturbing is that the full- 
sky map making algorithm is inserting significant extra 
large angle power into precisely those portions of the sky 
where we have the least reliable information. We are sim- 
ilarly concerned that the "galactic bias correction" of the 
ILC123, which inserts extra low-£ power into the region 
of the galactic cut, does so on the basis of the expecta- 
tions of a model whose assumptions (statistical isotropy 
and the particular ACDM power spectrum) are not borne 
out even after all the corrections are made. But even the 
first-year ILC, which did not have this bias correction, 
suffers from a surfeit of power inside the cut. 

Finally, although the maximum likelihood estimates 
of the low-i" Ce (actually C e ) quoted by WMAP in their 
three-year release appear to be in better agreement with 
the ACDM model than are the values one derives from 
cut-sky individual bands, or the cut sky ILC123, or the 
values quoted in the one-year WMAP release, neverthe- 
less, the angular correlation function that the new quoted 
Ci imply (as shown in Figure 5) appears to be in even 
worse agreement with the ACDM prediction than all the 
others. 



IV. CONCLUSIONS 

We have shown that the ILC123 map, a full sky map 
derived from the first three years of WMAP data like its 
predecessors the ILC1, TOH1 and LILC1 maps, derived 
from the first year of WMAP data, shows statistically sig- 
nificant deviations from the expected Gaussian-random, 
statistically isotropic sky with a generic inflationary spec- 
trum of perturbations. In particular: there is a dramatic 
lack of angular correlations at angles greater than sixty 
degrees; the octopole is quite planar with the three oc- 
topole planes aligning with the quadrupole plane; these 
planes are perpendicular to the ecliptic plane (albeit at 
reduced significance than in the first-year full-sky maps); 
the ecliptic plane neatly separates two extrema of the 
combined I = 2 and I — 3 map, with the strongest ex- 
trema to the south of the ecliptic and the weaker ex- 



trema to the north. The probability that each of these 
would happen by chance are 0.03% (quoting the cut-sky 
ILC123 S 1/2 probability), 0.4%, 10%, and < 5%. As 
they are all independent and all involve primarily the 
quadrupole and octopole, they represent a ~ 10 -8 prob- 
ability chance "fluke" in the two largest scale modes. To 
quote [7]: We find it hard to believe that these correla- 
tions are just statistical fluctuations around standard in- 
flationary cosmology's prediction of statistically isotropic 
Gaussian random ai m [with a nearly scale-free primordial 
spectrum]. 

What arc the consequences and possible explanations 
of these correlations? There are several options — they 
are statistical flukes, they are cosmological in origin, they 
are due to improper subtraction of known foregrounds, 
they are due to a previously unexpected foreground, or 
they are due to WMAP systematics. As remarked above 
it is difficult for us to accept the occurrence of a 10~ 8 
unlikely event as a scientific explanation. 

A cosmological mechanism could possibly explain the 
weakness of large angle correlations, and the alignment 
of the quadrupole and octopole to each other. A cosmo- 
logical explanation must ignore the observed correlations 
to the solar system, as there is no chance that the uni- 
verse knows about the orientation of the solar system nor 
vice-versa. These latter correlations are unlikely at the 
level of less than 1 in 200 (plus an additional indepen- 
dent w 1/10 unlikely correlation with the dipole which 
we have ignored). This possibility seems to us contrived 
and suggests to us that explanations which do not ac- 
count for the connection to solar system geometry should 
be viewed with considerable skepticism. 

In [16], we showed that the known Galactic foregrounds 
possess a multipole vector structure wholly dissimilar to 
those of the observed quadrupole and octopole. This ar- 
gues strongly against any explanation of the observed 
quadrupole and octopole in terms of these known Galac- 
tic foregrounds. 

A number of authors have attempted to explain the 
observed quadrupole-octopole correlations in terms of a 
new foreground [51-56]. (Some of these also attempted 
to explain the absence of large angle correlations, for 
which there are also other proffered explanations [57- 
63].) Only one of the proposals ([53]) can possibly ex- 
plain the ecliptic correlations, as all the others are extra- 
galactic. Some do claim to explain the less-significant 
dipole correlations. Difficulties with individual mecha- 
nisms have been discussed by several authors [56, 64-67] 
(sometimes before the corresponding proposal). Unfor- 
tunately, in each and every case, among possible other 
deficiencies, the pattern of fluctuations proposed is incon- 
sistent with the one observed on the sky. As remarked 
above, the quadrupole of the sky is nearly pure Y22 in the 
frame where the z-axis is parallel to u^ 2 ' 1 ' 2 ) (or any nearly 
equivalent direction), while the octopole is dominantly 
Y33 in the same frame. Mechanisms which produce an 
alteration of the microwave signal from a relatively small 
patch of sky — and all of the above proposals fall into this 
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class — are most likely to produce aligned Y20 and Y30. 
(This is because if there is only one preferred direction, 
then the multipole vectors of the affected multipoles will 
all be parallel to each other, leading to a Yg .) The au- 
thors of [55] manage to ameliorate the situation slightly 
by constructing a distorted patch, leading to an under- 
powered F33, but still a pure Y" 2 o- The second short- 
coming of all explanations where contaminating effect is 
effectively added on top of intrinsic CMB temperature is 
that chance cancellation is typically required to produce 
the low power at large scales, or else the intrinsic CMB 
happens to have even less power than what we observe. 
Likelihood therefore disfavors all additive explanations 
[56] (unless the explanation helps significantly with some 
aspect of structure seen at higher £). 

Explaining the observed correlations in terms of fore- 
grounds is difficult. The combined quadrupole and oc- 
topole map suggests a foreground source which form a 
plane perpendicular to the ecliptic, ft is clear neither how 
to form such a plane, nor how it could have escaped de- 
tection by other means. This planar configuration means 
that single anomalous hot or cold spots do not provide 
an adequate explanation for the observed effects. 

The final possibility is that systematic effects remain 
in the analysis of the WMAP data. Two systematic cor- 
rections have already been introduced between the one- 
year and the three-year data releases: the improved ra- 
diometer gain model, which includes effects correlated 
with the ecliptic and with the seasons (and hence the 
equinoxes); and the galaxy bias correction to the ILC. 
We have discussed the roles of both of these in reducing 
somewhat the statistical significance of the orthogonal- 
ity of the quadrupole and octopolc planes to the ecliptic. 
But other systematic effects may remain. For example, a 
suggestion has been made [8] that observed north-south 
power asymmetries [9] could be caused by a small error 
(within the la error bars) in the solar velocity with re- 
spect to the CMB. However, this would not explain the 
correlations in full sky maps noted in this paper [68]. 
Other possibilities have been put forward, for example a 
correlation between the scanning pattern and the beam 
asymmetry [69], but have not been carefully evaluated. 

If indeed the observed £ = 2 and 3 CMB fluctuations 
are not cosmological, there are important consequences. 
Certainly, one must reconsider [7] all CMB results that 
rely on low is, such as the normalization, A, of the pri- 
mordial fluctuations and any constraint on the running 
dn s /dlogk of the spectral index of scalar perturbations 
(which, as noted in [70], depended in WMAP1 on the 
absence of \ow-£ TT power). Moreover, the CMB-galaxy 
cross-correlation, which has been used to provide evi- 
dence for the Integrated Sachs- Wolfe effect and hence the 
existence of dark energy, also gets contributions from the 
lowest multipoles (though the main contribution comes 
from slightly smaller scales, £ ~ 10). Finally, note that, 
even though we have discussed only alignments at £ — 2 
and 3, it is possible - even likely - that the underly- 
ing physical mechanism does not cut off abruptly at the 



octopole, but rather affects £ = 4, 5, . . .. Indeed, sev- 
eral pieces of evidence, though less statistically convinc- 
ing than those discussed here, have been presented for 
anomalies at I > 3. If the values of Ce for £ > 3 are 
called into doubt, then so are the values of further cos- 
mological parameters, including the optical depth to the 
last scattering surface r, the inferred redshift of reioniza- 
tion, and the rms mass fluctuation amplitude cr 8 . 

Of even more fundamental long-term importance to 
cosmology, a non-cosmological origin for the currently- 
observed \ow-£ microwave background fluctuations is 
likely to imply further-reduced correlation at large an- 
gles in the true CMB. As shown in Section 3, angular 
correlations are already suppressed compared to ACDM 
at scales greater than 60 degrees at between 99.85% and 
99.97%C.L. (with the latter value being the one appropri- 
ate to the cut sky ILC123). This result is more significant 
in the year 123 data than in the year 1 data. The less 
correlation there is at large angles, the poorer the agree- 
ment of the observations with the predictions of generic 
inflation. This implies, with increasing confidence, that 
either we must adopt an even more contrived model of 
inflation, or seek other explanations for at least some of 
our cosmological conundrums. Moreover, any analysis 
of the likelihood of the observed "low-^ anomaly" that 
relies only on the (low) value of C2 (especially the MLE- 
inferred) should be questioned. According to inflation 
C2, C3 and C4 should be independent variables, but the 
vanishing of C{9) at large angles suggests that the differ- 
ent \ow-£ Ce are not independent. 

Another striking fact seen in Table VII is that the 
quadrupole and octopole of cut-sky maps, as well as 
S1/2 statistics, are significantly lower than in the full- 
sky ILC maps, both for year 1 and year 123. This seems 
to be because of a combination of effects — the "galac- 
tic bias correction" and the time-dependent radiometer 
gain model certainly, but perhaps also just the minimal 
variance procedure used to build the ILC. This reflects 
the fact that the synthesized full sky maps show corre- 
lations at large angles that are simply not found in the 
underlying band maps. This does not seem reasonable 
to us — that one starts with data that has very low cor- 
relations at large angles, synthesizes that data, corrects 
for systematics and foregrounds and then concludes that 
the underlying cosmological data is much more correlated 
than the observations — in other words that there is a 
conspiracy of systematics and foreground to cancel the 
true cosmological correlations. 

This strongly suggests to us that there remain serious 
issues relating to the failure of statistical isotropy that are 
permeating the map making, as well as the extraction of 
low-£ C e . 

At the moment it is difficult to construct a single co- 
herent narrative of the low I microwave background ob- 
servations. What is clear is that, despite the work that 
remains to be done understanding the origin of the ob- 
served statistically anisotropic microwave fluctuations, 
there are problems looming at large angles for standard 
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inflationary cosmology. 
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